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We consider Deeply Virtual Compton Scattering (DVCS) on nucleons and nuclei in the framework 
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of generalized vector meson dominance (GVMD) model. We demonstrate that the GVMD model 
provides a good description of the HERA data on the dependence of the proton DVCS cross 
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I. INTRODUCTION 



During the last decade, one of main focuses of hadronic physics has been the study of the 
hadronic structure using hard exclusive reactions, such as Deeply Virtual Compton Scat- 
tering (DVCS), 7*iV — > jN, and hard exclusive meson production (HEMP), 7*iV — > MN. 



1 Sy£,b-S e £ the subjec ; of intensive theoretical and experimental investiga- 



te] 
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In addition, there were investigated "inverse" 



Thesepro cesses 

tions 0,y,y,y, 

hard exclusive reactions such as 'yN — > 7*iV — > l + l~N 
and "w-channel" reactions such as 7*7 — > tttt 

The interest to the DVCS and HEMP reactions is motivated bv the fact that in the 



and ttN 7 *iV l+l~N [16|, 



?a Dy 

m 



Bjorken limit (large Q 2 ), the corresponding amplitudes factorize [7|, [8| in convolution of per- 
turbative (hard) coefficient functions with nonperturbative (soft) matrix elements, which 
are parameterized in terms of generalized parton distributions (GPDs). GPDs are univer- 
sal (process-independent) functions that contain information on parton distributions, form 
factors and correlations in hadrons. GPDs also parameterize parton correlations in matrix 
elements describing transitions between two different hadrons, which appear in reactions 
such as e.g. 7*p — > 7r + n Q, Q- 

While the description of DVCS and HEMP based on the factorization approach is most 



general, in experiments the values of the virtua 
the validity of the factorization theorem 0, E 



ities Q 2 are 



DC 
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ow the range required for 



251 ] . Hence, contributions of 



higher- twist effects might be substantial (it is an open issue how large these effects are), 
which will affect the extraction of GPDs from the data. Therefore, it is important to have 
an effective model for the DVCS and HEMP amplitudes, which would interpolate between 
the photoproduction (Q 2 w 0) and deep inelastic (Q 2 ~ O(10) GeV 2 ) regimes. 



In this paper, using the generalized vector meson dominance (GVMD) model [26l. l27l. |28|. 
which is consistent with perturbative QCD at small transverse distances j^, we derive 
expressions for the amplitudes of DVCS on nucleons and nuclei, which are valid at high 
energies and which are applicable over a wide range of Q 2 . We show that the resulting cross 
section of DVCS on nucleons compares well to the HERA data Q Q- In particular, the 
dependence of the DVCS cross section on Q 2 , W (at Q 2 =4 GeV 2 ) and t are reproduced rather 
well; the W-dependence of the cross section at Q 2 =S GeV 2 is somewhat underestimated, 
which can be interpreted as due to the onset of the hard regime beyond the soft dynamics 
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of the GVMD model. 

We also estimate the relative contribution of 
higher-twist corrections in perturbative QCD 
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/(^-corrections, which correspond to the 
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34|. We show that these corrections 



are large: the contribution of the l/Q 2 -corrections to the DVCS amplitude at t — t min is 
20% at Q 2 = 2 GeV 2 , 11% at Q 2 = 4 GeV 2 and 6% at Q 2 = 8 GeV 2 ; the contribution of 
the l/Q 2 -corrections to the t- integrated DVCS cross section is 56% at Q 2 = 2 GeV 2 , 32% 
at Q 2 = 4 GeV 2 and 17% at Q 2 = 8 GeV 2 . 

We also make predictions for the DVCS cross section on nuclear targets, which are relevant 
for the physics program of the future Electron-Ion Collider. We predict significant nuclear 
shado wing , which matches well predictions of the leading-twist nuclear shadowing in DIS on 



35] 



nuclei 

The hypothesis of vector meson dominance (VMD) assumes a definite relation be- 
tween the amplitude of the photon (real or virtual)-hadron interaction, A(j* T + T — > . . . ), 
and a linear combination of the amplitudes of the corresponding strong production by trans- 
versely polarized vector mesons, A(V tr + T — > . . . ), 

A(jt + T^...)= £ A(V tT + T^...), (1) 

where fv is the coupling constant determined from the V — > e + e~ decay; my is the vector 
meson mass; Q 2 is the virtuality of the photon; T denotes any hadronic target. Note that 
Eq. ([I]) is written for the transversely polarized photons. In Eq. ([T|), we took into account 
only the contribution of the p°, u and mesons. 

The VMD model and its generalizations explain a large wealth of data on the real and 
virtual (Q 2 < 1 GeV 2 ) photon-hadron scattering, which include the pion electric form fac- 
tor, total cross sections of photon-nucleon and photon-nucleus scattering (inclusive structure 
functions), exclusive production of vector mesons on nucleons and nuclei, exclusive produc- 



tion of pseudoscalar mesons, for a review, see [371 ]. 



As the virtuality of the photon increases, Q 2 > 1 GeV 2 , the simple VMD model, see 
Eq. (CQ), becomes inadequate since it leads to the violation of the approximate Bjorken 
scaling. In order to restore the approximate Bjorken scaling, the simple VMD model can be 



generalized 38]. This can be done using the model- independent method of mass- dispersion 



representation for the virtual photon-hadron scattering amplitude 39]. 



In order to illustrate the approach, let us consider the forward virtual photon-hadron 



scattering amplitude. The dispersion representation for the imaginary part of A(^ r + T — > 
7 t * r + T) reads 

r, ,/ * m * ^ /" dM 2 M 2 dM' 2 M' 2 e e 

where cr^v" is the V + T — ► V' + T scattering cross section (spectral function) which weakly 
depends on the masses M and M'. The main idea of the GVMD model 23, Um, [28( is to 
approximate Eq. (j2J) by an infinite series of (ficticious) vector mesons of ever increasing mass, 
allowing for both diagonal (V + T — ► V + T) and non-diagonal (V + T — > V' + T) transitions. 
The role of the non-diagonal transitions is to partially cancel the diagonal transitions so that, 
effectively, uvv 1/M 2 for large M 2 . This softens the spectral function and leads to the 
approximate Bjorken scaling of the inclusive structure function F 2 (x,Q 2 ), see [29!] for the 
discussion. 

In the language of the color dipole model, the fact that oyy oc 1/M 2 for large M 2 
means that besides dipoles of large transverse sizes, the virtual photon also contains small 
transverse-size dipoles. The latter fact is called color transparency. 

One should note that, while the simple vector meson dominance model fails to reproduce 
the approximate scaling of the inclusive structure function F 2 (x, Q 2 ) (see above), the simple 
VMD model predicts the correct Q 2 -behavior of cross sections of exclusive reactions, such 



as 7*p — > n^n 4U|, |4l|]. This is also true for DVCS: Even the simple VMD model provides 
the correct Q 2 -behavior of the DVCS cross section (up to logarithmic corrections). 

The structure of this paper is the following. In Sect. [Til we explain main assumptions of 
the GVMD model, which we further generalize to take into account a non-zero momentum 
transfer t ^ 0. We derive the expression for the amplitude of DVCS on the nucleon and make 
predictions for the DVCS cross section. We demonstrate that the GVMD model provides a 
good description of the HERA data on the W, Q 2 and t-dependence of the cross section of 
DVCS on the proton 30, 3l|. In this section, we also estimate l/Q 2 -corrections to the DVCS 



amplitude and the DVCS cross section. Predictions for the nuclear DVCS amplitude and 
for the cross section of DVCS on nuclei in the collider kinematics are presented in Sect. IIHI 
In Sect. ITVl we summarize and discuss our results. 
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FIG. 1: The DVCS amplitude in the GVMD model, see Eq. ©. 



II. DVCS ON THE NUCLEON 



27 



n this section, we extend the generalized vector meson dominance (GVMD) model 26, 



28] to the off-forward case and apply it to Deeply Virtual Compton Scattering (DVCS) 



on the nucleon. 



A. DVCS amplitude 

The GVMD model assumes that the virtual (real) photon interacts with the hadronic 
target by fluctuating into a coherent and infinite sum of ficticious vector mesons V n . Then, 
the DVCS amplitude at the photon level, A(j*P ~~ > IP), can De graphically presented as 
depicted in Fig. CD 

In the GVMD model, the DVCS amplitude for transversely polarized virtual photons 
reads 

•A(t£ + p-7 + p)= E 7m^^'^V- ( 3 ) 

n,n'=0 Jn 1Un + ^ Jn ' 

where W 2 = (p 7 * + p) 2 with p 7 * the momentum of the initial photon and p the momentum 
of the initial proton; t = (p' — p) 2 with p' the momentum of the final proton. The masses 
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M n and the coupling constants /„ are connected by the following relations, 

M 2 f 2 

wrfr {l + 2n) ' (4) 

where M = m p and / — fp refer to the physical p° meson. 

Note that the relation of the vector mesons V n conventionally used in the GVMD model 
to physical J p = 1~ vector mesons found in the Review of particle physics [4^] is not direct. 
The motivation for this is that while the vector meson masses are known with a reasonable 
accuracy up to M p > 2 GeV, there is no accurate data on the partial decay width T e + e - 
for mesons heavier than p(1450). On the other hand, the parameterization (J4j) provides 
reasonable results for physical observables. One can check that the linear n-dependence of 
the ratio M 2 /Mq in Eq. (jlj) is confirmed experimentally for large-n 42j. However, the slope 
of the n-dependence is underestimated by approximately a factor of two. 

It is important to point out that, at high energies, the DVCS cross section at the photon 
level is dominated by the contribution of the transversely polarized virtual photons due to 
the helicity conservation 3l| . In the language of the color dipole model, this dominance is 
explained by the dominance of the large transverse-size dipoles over the small transverse-size 



dipoles, see 



131 ] for the discussion. Therefore, Eq. ([3]) gives the complete description of the 



DVCS amplitude. 

The quantity S n n /(t) is the V n>tT + p — > V n / )tr + p scattering amplitude, see Fig. [IJ The 
matrix E ra „/(t) is assumed to have a tri-diagonal form with the following non-zero elements, 

Z n ,n(W,t) = ia pp (W 2 )(l-ir l )F n (t), 

£ n>n+1 (w;t) = x n+hn (w,t) = -~- AL 



I -2^)E BfB (W;t) 



2M n+1 V ~~M 2 ) -«."v"'-/' ( 5 ) 
where <7 PP is the p meson-proton cross section for the transversely polarized meson; 77 is the 
ratio of the real to imaginary parts of the p meson-proton scattering amplitude; 5 = 0.2 is 
the parameter of the model. 

The function F n (t) models the t-dependence of E njn /(t), which goes beyond the original 



formulation of the GVMD model 



26 



27 



, which addressed only the forward t — limit. 



In our analysis, we use the following form of F n (t), 



F n (t) = exp 



n 



-B,+ 



n 



n+l 



-Bo 



(6) 



where B\ = 11 GeV and B2 = 4.3 GeV . The choice of the slopes B\ and B2 is motivated 



as follows. 
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For the moment, let us replace the final real photon by the p meson. In the photopro- 
duction limit, the 7p — > pp cross section measured at HERA by the HI collaboration was 
fitted to the form exp(-B\t\) with the slope B = (10.9 ± 2.4 ± 1.1) GeV" 2 43]. The ZEUS 



44J. 



measurement gives essentially the same value of B 

In electroproduction, the slope of the exponential fit to the ^*p — > pp cross section is 
much smaller than in photoproduction: It decreases from B = (8.0 ± 0.5 ± 0.6) GeV -2 at 



Q 2 = 1.8 GeV 2 to B = (4.7 ± 1.0 ± 0.7) GeV" 2 at Q 2 = 21.2 GeV 2 [45j. 

This decrease of the slope of the t-dependence with increasing Q 2 is effectively parame- 
terized by Eq. ([6]) as a decrease of the slope with the increasing number of the vector meson 
n. Indeed, close to the photoproduction limit, the dominant contribution to the sum in 
Eq. ([3]) comes from the n = term. In the opposite limit of large Q 2 , terms with large n, 
up to M 2 ~ Q 2 , contribute to the sum. Choosing Q 2 = and Q 2 = 21.2 GeV 2 as reference 
points, we determine the values of the slopes Bi and B 2 in Eq. ([6]). 

It is clear from the above discussion that while the value of the slope B\ in Eq. ([6]) is 
model-independent, the value of the slope B 2 is somewhat more uncertain. We have chosen 
not to introduce Q 2 -dependent slopes B\ and B 2 since this would contradict the spirit of 
the VMD model: The W and t-dependence of the DVCS amplitude is determined solely by 
the vector meson-proton scattering amplitudes; the vector meson propagators provide the 
Q 2 -dependence. 

It is important to note that the non-diagonal terms, £„,„' with n ^ n', are essential in 
the GVMD model: The infinite series in Eq. ([3]) would have been divergent without the 
non-diagonal transitions. Also, the non-diagonal terms provide the correct scaling of the 
total 7*p cross section. 

Using Eqs. (jl]) and (Q, the DVCS amplitude in Eq. ([3]) can be written in the following 
form 



F n (t) 



x 



^ "" / ^(^ + l + 2n)(^+3 + 2n) 

Q 2 (, ,r 1 + U 



1 + 2m 2 (l + 5) (3 + 2n) \ + 45 (1 + 2n) 



(7) 



Equation (J7|) involves four quantities, f p , a pp , rj and 5, which are known with a certain degree 
of uncertainty. One can reduce this uncertainty by expressing the DVCS amplitude in terms 
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of the total ■yp cross section, 



J o 



pp 



n=0 



+ 2n)(3 + 2n) 



T 2 

J o 



o PV {W 2 



(8) 



Therefore, the final expression for the DVCS amplitude reads 

oo 

Q 2 



F n (t) 



X 



1 + 



+ l + 2n)(^ + 3 + 2n) 
1 + n 



1 + 45- 



(9) 



2m 2 p (l + 5)(3 + 2n) \ (l + 2n)/. ' 

One should also note that another advantage of expressing the DVCS amplitude in terms 
of a2ot i s that Eq. ([9]) effectively takes into account the contributions of the u and <fr vector 
mesons, which enter through the phenomenological parameterization of a^ t . 



In our analysis, we use the ZEUS parameterization of cr^fW 2 ) [47]] 



57W°- 2 + 121W 



-0.716 



(10) 



where the cross section is in fib and W is in GeV. 



The ratio of the real to imaginary parts o: 



is found using the Gribov-Migdal relation 



48 



the V + p —> V + p scattering amplitude, rj, 



rr p 

-- « 0.16 
2 2 



(11) 



where p = 0.2 was used, which corresponds to the power of the H^-dependence of cr^ t {W 2 ) 
at large W in Eq. (JTOJ). 



26 



28 



46] 



The remaining parameter in Eq. Qj is 6, for which we use 5 = 0.2 
However, the exact numerical value of 5 affects weakly our numerical predictions. 

One of simplest DVCS observables is the skewing factor R, which is defined as the ratio of 
the DVCS to the DIS amplitudes and which was recently extracted from the HERA DVCS 
and DIS data [49(, 



QmAdfr + p — > 7 + p) 



I* 



SmA(yt T +p-> 7tr + P)\t=o 



(12) 



Note that we generalized the ratio R originally defined at t = t m ; n 49] to any value of t. At 
high energies, the minimal momentum transfer |t m i n | ~ x 2 B m 2 N w 0, where Xb is the Bjorken 
variable; is the nucleon mass. 
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FIG. 2: The GVMD prediction for the ratio R(t) of the DVCS and DIS amplitudes, see Eq. (fT2|) . 
as a function of Q 2 for three values of t. 



The GVMD model makes an unambiguous prediction for the ratio R, 

^ + 1 + 2n)(S + 3 + 2n) L 2m 2 (3 + 2n) V (1 + 2n)J 



E (S+1+ ^| +3 ,^ ^H- 



•p 




Figure [2] presents the GVMD predictions for the ratio R as a function of Q 2 for three 
values of t: t = t min » (evaluated with x B = 0.001), t = -0.1 GeV 2 and t = -0.2 GeV 2 . 
Note that in the GVMD model, the ratio R does not depend on W or xb at given Q 2 and 
t. 

A comparison of the solid curve in Fig. [2] to the experimental results for the ratio R, see 
Fig. 4 of Ref. j49j], reveals that the GVMD model provides a good description of the data for 
Q 2 < 5 GeV 2 . For higher values of Q 2 , the GVMD model overestimates the experimental R. 
Therefore, the GVMD model and similar models can be used to reliably determine DVCS 
observables and generalized parton distributions at Q 2 of the order of a few GeV 2 . This can 
be used as an input for QCD evolution to higher Q 2 scales. An example of such an approach, 
which uses the align-jet model to construct input GPDs and which excellently compares to 
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the HERA data on the DVCS cross section and on the ratio R, was worked out in [5 



B. DVCS cross section 

The DVCS amplitude in Eq. (Q is normalized such that in the Q 2 — > limit, the imagi- 
nary part of the 7p — > ■yp amplitude is equal to the total photoabsorption cross section, see 
Eq. (jSJ). With such a normalization, the differential and integrated DVCS cross sections at 
the photon level read 

dt l07T 

^DVCS(^Q 2 ) = T7- / dtlA^r+P^l+PW- (14) 

lO 71 " J -I GeV 2 

where £ m j n ~ —x 2 B m 2 N (t m i n ~ in the HERA kinematics). 

In order to compare the GVMD model predictions to the data on the DVCS cross section 



at the photon level 30, |31[, one needs to make sure that one compares the same quantities. 
Using the classic result of L. N. Hand [51J], one can readily see that the HERA DVCS cross 
section at the photon level is indeed a properly defined and normalized cross section of the 
7*p — > 7p reaction. 

As a byproduct of the above mentioned exercise, one establishes the connection between 
the GVMD and GPD descriptions of the DVCS cross section: 

i^(7;+p-7+p)i 2 = 4 e !^-_ ((1 - e) m\ 2 + \n\ 2 ) - e(n*£ + us* 



+ H . £+w .)_i m _ + a-^i£pj, us) 

where %b is the Bjorken variable; £ = £b/(2 — xs); e 2 = 4x 2 3 m 2 v /Q 2 . The quantities 
7i and ^ are the so-called Compton form factors of the corresponding proton GPDs 12]. 
It is important to have the connection between the GVMD-based and the GPD-based de- 
scriptions of the DVCS cross section since the both approaches have an overlapping region 
of applicability, namely, 1 < Q 2 < 5 GeV 2 . 

The simple expression for the DVCS amplitude in the GVMD model ([9]) allows one 
to examine the relative contribution of l/Q 2 -corrections, which correspond to higher- twist 
corrections in perturbative QCD. To this end, let us expand the DVCS amplitude in Eq. ([9]) in 
terms of l/Q 2 and let us call the leading contribution, which behaves as l/Q 2 , A^ij^+p — > 
7 + p). The corresponding t- integrated cross section is denoted as crpy CS . 
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TABLE I: The l/(5 2 -corrections to the DVCS amplitude and to the i-integrated DVCS cross section 
as functions of Q 2 , see Eq. (fl"6j) . 



Q 2 [GeV 2 ] 




R? T (Q 2 ) 


2 


0.20 


0.56 


4 


0.11 


0.32 


8 


0.058 


0.17 



We quantify the contribution of l/Q 2 -corrections to the DVCS amplitude and to the 
DVCS cross section by introducting the ratios -R^mpi an d R^ T , 

pHT (r) 2, _ i _ -4 LQ (7 t * r +P^7 + P)|^ min 

i?HT (Q2) = x _ ^D§CS (16) 
0"DVCS 

The ratios -R^mpi an d -^a T as functions of Q 2 are summarized in Table fl] Note that these 
ratios do not depend on W in the chosen model. 

As one can see from Table [J the 1 /Q 2 -corrections are large. Moreover, R^ T (Q 2 ) > 
2R^ pl (Q 2 ) due to the enhancement of the heavy vector meson contributions to ctdvcs be- 
cause of the decreasing slope of the t-dependence of the DVCS amplitude with increasing n, 
see Eq. ©. 



C. Comparison to the HERA DVCS data 



Using Eqs. ([HD and (FH|) . we make predictions for the DVCS cross section and compare 



our findings to the HERA data 



30 



3l|. 



Figure [3] presents the TV-dependence of the DVCS cross section at Q 2 = 4 GeV 2 and 
Q 2 = 8 GeV 2 . The solid curves correspond to the GVMD calculations; the experimental 



points are those from the HI 



3JJ and ZEUS 30j measurements. The error bars correspond 



to the statistical and systematic errors added in quadrature. The ZEUS data taken at 
Q 2 = 9.6 GeV 2 have been interpolated to Q 2 = 8.0 GeV 2 using the fit to the Q 2 -dependence 
of ctdvcs, cloves ~ l/(Q 2 T with u = 1.54 y. 

We shall discuss the left and right panels of Fig. [3] separately. As seen from the left panel 
of Fig. [3l the GVMD model reproduces both the absolute value and the ly-dependence of 
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20 40 60 80 100 120 140 160 20 40 60 80 100 120 140 160 



W [GeV] W [GeV] 

FIG. 3: The DVCS cross section as a function of W. The GVMD model results (solid curves) are 
compared to the HI [3j| and ZEUS [30| data. The error bars correspond to the statistical and 
systematic errors added in quadrature. 



ctdvcs sufficiently well. The latter fact signifies that, at Q 2 = 4 GeV 2 , the DVCS cross 
section is still dominated by soft physics. At Q 2 = 4 GeV 2 , the VU-behavior of ctdvcs is 
consistent with that predicted by the GVMD model, ctdvcs ~ W 0A . 

Turning to the right panel of Fig. [3], we observe that while the GVMD model compares 
fairly with the HI data, the model underestimates the slope of the VF-dependence of ctdvcs 
for the ZEUS data set, which has smaller error bars. In particular, the predicted odvcs ~ 
W 0A behavior is much slower than that given by the fit to the ZEUS data points, ctdvcs ~ W s 



with 5 = 0.75 ± 0.15 30j. This indicates the onset of the hard regime in the total DVCS 
cross section at Q 2 = 8 GeV 2 , where the GVMD model becomes inadequate. 

Figure H] presents the Q 2 -dependence of the DVCS cross section at W = 82 GeV. The solid 
curve corresponds to the GVMD model; the experimental points come from the HI [31] and 



ZEUS 



30j experiments. The error bars correspond to the statistical and systematic errors 



added in quadrature. The ZEUS data taken at W = 89 GeV have been extrapolated to 



W = 82 GeV using the fitted IU-dependence of a DV cs, odvcs ~ IU a75 [301 ]. 

One sees from Fig. 0] that the GVMD model reproduces the Q 2 -dependence of aovcs 
over a very wide range of Q 2 , 3 < Q 2 < 85 GeV 2 . This is quite a remarkable result that 
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FIG. 4: The DVCS cross section as a function of Q 2 . The GVMD model result (solid curve) is 



compared to the HI I3J and ZEUS 
systematic errors added in quadrature. 



301 ] data. The error bars correspond to the statistical and 



the model, which was initially developed for photoproduction and was later extended to 
electroproduction with Q 2 of the order of a few GeV 2 , provides a quantitative description for 
such large values of Q 2 . In other words, at fixed W, the GVMD model correctly reproduces 
the Q 2 -scaling of ctdvcs- 

Figure [5] presents the comparison of the t-dependence of the GVMD model ca 
(solid curves) to the HI data on the differential DVCS cross section derives jdt 
error bars are the statistical and systematic errors added in quadrature. 

As one sees from Fig. the GVMD model describes danves/dt well. This result is not 
quite trivial. In order to achieve this within the framework of the GVMD model, one has to 
assume that either all vector meson V^-nucleon cross sections have the same Q 2 -dependent 
slope of the ^-dependence or that the slope decreases with increasing n, see Eq. ([6]). While 
the value of the slope B\ — 11 GeV -2 is fixed by photoproduction of p mesons, the value of 
the slope B 2 is model-dependent. The values B 2 — 4 5 GeV~ 2 provide a good agreement 
with the HI data (see Fig. E]), which were fitted to the exponential form, da^yes/dt ~ e^- 8 '*' 
with B = (6.66 ± 0.54 ± 0.43) GeV" 2 at Q 2 = 4 GeV 2 and B = (5.82 ± 0.59 ± 0.60) GeV~ 2 
at Q 2 = 8 GeV 2 H. 
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FIG. 5: The differential DVCS cross section daj)vcs/dt as a function of t. The GVMD model 



calculation (solid curves) is compared to the HI data 
systematic errors added in quadrature. 



311 ] . The error bars are the statistical and 



D. DVCS cross section in Jefferson Lab kinematics 



The ep —>■ epj cross section in the DVCS regime was recently measured by the Hall A 
collaboration at Jefferson Laboratory (JLab) 24j. The cross section involves the contri- 
butions of the Bethe-Heitler (BH) amplitude squared, the DVCS amplitude squared and 
the interference of the BH and DVCS amplitudes. Based on the kinematics of the exper- 
iment, in the analysis of the data, the contribution of the DVCS amplitude squared was 



12], which allowed for the extraction of 



neglected compared to the other two contributions 
the so-called Compton form factors of the proton. 

In this subsection, we check the validity of the assumption that the contribution of the 
DVCS amplitude squared is negligibly small by explicitly calculating the DVCS cross section 
within the GVMD model in the Jefferson Lab kinematics. The DVCS cross section at the 
lepton level reads, see e.g. [121 ]. 

[l-y + y 2 /2) 1 da mC s(W,Q 2 ,t) 



d 4 a 



a:, 



(17) 



dQ 2 dxBdtd(f) itQ 2 xb 27T dt 

where a e . m . is the fine-structure constant; is the angle between the lepton and production 
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planes; d<TDVcs(W, Q 2 , t)/dt is the DVCS cross section at the proton level defined by Eq. (JHj). 
The extra factor 1 / (2ir) in the right-hand side of Eq. (TPTj) takes into account the fact that 
the integration over the angle is included in the definition of cZcr DVCS (TV, Q 2 , t)/dt. The 
DVCS cross section does not depend on <fr when one neglects the photon helicity changing 
transitions [111 ]. 



Using Eqs. (|9]) and (JHJ), we evaluate the DVCS cross section at the lepton level 
d 4 a/(dQ 2 dxBdtd(p) in the kinematics of the Hall A experiment, E = 5.75 GeV (the en- 
ergy of the lepton beam), Q 2 = 2.3 GeV, t = —0.28 GeV and x B = 0.36, 

' ° 0.0022 nb/GeV 4 . (18) 



dQ 2 dxsdtd(f) 

This value is an order of magnitude smaller than the sum of the BH and interference con- 
tributions to the ep — > epy cross section, which confirms the assumption that, in the JLab 
kinematics, the contribution of the DVCS amplitude squared to the ep epj cross section 
can be safely neglected. 

III. DVCS ON NUCLEI 

In this section, we derive the expression for the DVCS amplitude on a nucleus using 
the GVMD model for the photon-nucleon interactions and the generalized Glauber formal- 
ism [3_7j in order to account for the multiple rescattering of the vector mesons inside the 
nucleus. Using the obtained amplitude, we make predictions for the nuclear DVCS cross 
section at the photon level in the collider kinematics. 

At high energies, in the GVMD model, photons (real and virtual) interact with hadrons 
by fluctuating into an infinite sum of vector mesons. When the involved hadron is a nucleus, 
each vector meson undergoes multiple interactions with the nucleons of the nucleus, which 
leads to the attenuation (decrease) of the vector meson-nucleus cross section compared to 
the sum of free vector meson-nucleon cross section. As a consequence, the resulting photon- 
nucleus cross section is smaller than the sum of the corresponding photon-nucleon cross 
sections. This phenomenon is called nuclear shadowing. It has been observed in various 



37, \5j% for review. 



reactions with nuclei induced by real and virtual photons, see 

In the GVMD model, the nuclear DVCS amplitude can be organized as a multiple scat- 
tering series (Glauber series), where each term corresponds to the number of interactions of 
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. n m .' %n my 

% <. n' y n' n" / 





FIG. 6: A schematic representation of the multiple scattering (Glauber) series for the nuclear 
DVCS amplitude in the GVMD model. The dashed lines correspond to vector mesons; the solid 
lines correspond to nucleons; the ovals with legs correspond to the final and initial nucleus. 



the vector mesons with the nucleons. This is schematically presented in Fig. [6j where the 
interactions with one, two and three nucleons are depicted. The dashed lines correspond 
to the vector mesons (note that the GVMD model allows for non-diagonal vector meson- 
nucleon transitions); the solid lines correspond to the nucleons involved in the interactions; 
the nuclear part is denoted by ovals with legs corresponding to the initial and final nucleus. 



Using the standard technique 37|, one can readily write down the expression for the 
nuclear DVCS amplitude A(Y + A -> 7 + A), 



A(y + A^ 1 + A) = J2 



Ml 



n,m=0 

A(A - 1 



AF A {t)E n>m {W,t) 



poo 

dzi / dz 2 I d 2 be icit - h p{b ) z l )p{b ) z 2 )e lZl ^'- kv - ) 



■2i 

( A — 2 ~ 
x £ nin /(W,Zi)(5 n> « —J: n/y >(W,z')e(z 1 <z'< z 2 )p(b,z') 



+ ... )t n ,, tm (W,z 2 )e 



2i 

iz 2 (k Vm -k 1 ) 



(19) 



where A is the number of nucleons in the nucleus (we do not distinguish protons and neu- 
trons); .Fa(£) is the nuclear form factor (-Fa(O) = 1); p(r) is the density of nucleons in the 
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nucleus [53]; b is the two-dimensional vector (impact parameter) in the plane perpendic- 
ular to the direction of the incoming photon, whose momentum is assumed to be along 
the ^-direction; Z{ are longitudinal positions of the nucleons of the nucleus involved in the 
interaction; q t is the transverse component of the momentum transfer. Since we neglected 
the ^-dependence of the elementary vector meson-nucleon scattering amplitudes compared 
to the steep t-dependence of the nuclear form factor, all scatterings of the vector mesons in 
Eq. (fT9|) occur at the same impact parameter b. 
In Eq. (DED, 

Km(W,z) = e iz ^- k ^X n>m (W,t = 0) , (20) 

where £„ jm is defined by Eq. (EJ). In Eqs. ( TT9l and (1201) . the exponential factors (except for 
the exp(iq t ■ b) factor) arise due to the non-zero longitudinal momentum transfer associated 
with non-diagonal in mass transitions. At high energies, 

k r ~ k Vn = V^+Ql - V" 2 ~Ml = Q2 + M " = x B m N fl + ^£ 

kv m ~k,= V^M! - v> = « = -XBm N ^ , 

M 2 - M 2 

kv n - k Vm = m ^ n , (21) 

where v is the energy of the incoming virtual photon in the laboratory reference frame; v' 
is the energy of the final real photon. We also used that v 1 — v + t/(2m^) w v for the small 
momentum transfer t. 

In Eq. (j!9p . the first term corresponds to the left graph in Fig. [61 which describes the 
interaction with one nucleon of the nucleus (the Born term). The second term in Eq. (IT5|) 
corresponds to the middle graph in Fig. which describes the interaction of hadronic 
fluctuations of the involved photons with two nucleons of the nucleus. Those nucleons are 
located at the points f\ = (b,z\) and r 2 = (b,z 2 ). This graph leads to the attenuation 
(nuclear shadowing) of the Born term. The third term corresponds to the interaction with 
three nucleons of the nucleus. The dots in Eq. (Tl9|) denote higher rescattering terms not 
shown in Fig. 

Equation <HM is rather general and, because of the non-diagonal V n V n > transitions, 
the direct calculation of the nuclear DVCS amplitude for heavy nuclei using Eq. (fT9l) is 
impossible. Therefore, for our numerical predictions, we make an approximation and ignore 
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the non-diagonal transitions for the interactions with three and more nucleons (this does 
affect the convergence of the series), 

A — 2 

S n ',n" %n',n»(W, z')Q(z 1 < Z < Z 2 )p(b, z') + . . . 



2i 



S n 'n" 1 



A-2 

2i 



E n >,n>(W,t = Q)e(z 1 <z'< z 2 )p(b,z') 



X , t/p -^ PP (W 2 )(l-i V )f:?dz>p(b,z>) 



(22) 



where in the last line we used the large- A approximation. Therefore, Eq. ( fl9l) now reads 

Ml 



n,m=0 

A(A-V r ' x 



f n M% + Q*f n 



AF A (t)E n , m (W,t) 



2i 



oo J Z\ 



d Zl I dz 2 I d 2 be i ^ b p(b,z 1 )p(b,z 2 )e tZl(h -' 



iz\ (k * -k Vn ) 



x £ n>n ,(W, Zl ) e ~^ {w2) ^ ^ dz ' p{b > z,) t n ,, m (W, z 2 )e iz ^- k ^ 



(23) 



For comparison of nuclear shadowing in DVCS and DIS, we also give the expression for 
the forward nuclear DIS amplitude, which can be readily obtained from Eq. ([231 . 



-4(7* + A - 1 * + A) lt=0 = 
A(A - 1 



e Ml e Ml 



2% 



n,m=0 
oo p oo 



AX n , m {W,t = 0) 



f n Ml + Qif m Ml + Q* 

POO p 

dzi / dz 2 / d 2 bp{b ) z l )p{b ) z 2 )e iz ^'- kv - ) 



— OO J 21 

A 



x t^W.z^e--^ 2 ^- 1 ^^ 



(24) 



We quantify predictions of the GVMD model for the nuclear DVCS and DIS amplitudes 
by considering the ratios R^ pl , R^ P \ and #£S>1> 

Im SmA(Y + A -> 7 + A) 



;,m|>1 ~ 3m^ Born (7* + A^7 + A) 
Ke^l(7* + A -> 7 + A) 



pRe 



DIS 

ampl o 



:e^t Born (7* + A -»• 7 + A) ' 
3m„4( 7 * + A -> 7* + A)i t=0 



where -4(7* + A -> 7 + A) is the nuclear DVCS amplitude of Eq. (T23J; ^l Born (7* + A -> 7 + A) 
is the first term (Born contribution) of Eq. ([231) ; ^l Born (7* + A — > 7* + A)| t= o is the first 
term of Eq. ([26]) . Note that -Rampi * s nothing but the the ratio of the nuclear to the nucleon 
inclusive structure functions, R^pi = F 2 a(x,Q 2 )/[AF 2N (x,Q 2 )}. 
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Figure [7] presents the ratios -R^ pl (solid curves), R^pi (dotted curves) and -R^pi (dot- 
dashed curves) at Q 2 = 1 GeV 2 and t = t min « —x 2 B m 2 N as functions of xb- The left panel 
corresponds to the nucleus of 40 Ca; the right panel corresponds to 208 Pb. 

Let us now discuss the results presented in Fig. [7J in detail. The solid and dot-dashed 
curves coincide for xb < 0.01 and deviate only slightly for 0.01 < xb < 0.1, which means 
that the amount of nuclear shadowing is the same in the imaginary parts of the DVCS and 
DIS amplitudes. This observation agrees with the results obtained within the framework of 
a different approach to nuclear GPDs at small-x^, when the latter are modeled using the 



align-jet model for the nucleon GPDs and a parameterization of usual nuclear PDFs 54]. 



Moreover, the amount of nuclear shadowing predicted by our calculations in the GVMD 
model matches very well the leading-twist predictions for F 2 a(x, Q 2 )/[AF 2 n(x, Q 2 )] made 
at somewhat higher Q 2 35J. This is a consequence of the fact the GVMD model predicts 
a significant amount of inclusive diffraction in 7*p scattering at all Q 2 , which controls the 



size of nuclear shadowing in the leading-twist theory of nuclear shadowing 35 1. 

As the value of xb is increased (at fixed Q 2 ), the shadowing correction decreases due to 
the decrease of a pp (W) and due to the increasingly destructive role of the e lZl ^*~ kVn ' and 
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FIG. 8: The ratio R^ pl at Q 2 = 1 GeV 2 as a function of t. The left panel is for 40 Ca; the right 
panel is for 208 Pb. 



e iz 2 (fcv m -fc-r) factors in Eq. 

For the ratio R^ p \ of the real parts (dotted curves), at small-x^, the shadowing correction 
is approximately two times as large as for the ratio of the imaginary parts because of the 
fact that (1 — rf) 2 = 1 — rf — 2irj, see Eqs. (E]) and (!23|) . As increases, the real part of 
the shadowing correction receives a large contribution from the e %z ^ k i*- k v n ) anc i e ^2(kv m -kj) 
factors, and, as a result, it steadily grows and becomes larger than the Born contribution. 



541 ] . Note, however, that since the 



This behavior of -R^Li ^ s similar to that observed in 
effect of t m ; n (the factor FAitmm) in the Born term) was not included in the analysis of [54j, 
the agreement could be coincidental. 

We also examined nuclear shadowing in DVCS with nuclear targets as a function of the 
momentum transfer t. Figure [8] presents -R^mpi at Q 2 = 1 GeV 2 as a function of t. Figured] 
demonstrates that the shadowing correction to the DVCS amplitude has the t-dependence 
which is slower than that of the Born term. As one increases \t\, the negative nuclear 
shadowing correction decreases slower than the Born term, which leads to a decrease of 
i?ampi- The ratio -R^ pl follows the similar trend. 

Equation (j23|) presents the 7*A — > scattering amplitude as a function of W, Q 2 and 
t. It also allows for the representation of the scattering amplitude as a function of W, Q 2 
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and b, where b is the impact parameter conjugate to q t 

Ml 



x « f 2 



n,m=0 Jn 1V1 n^^ /i 



/oo 
dze iXBmNZ p(b, z)E njfn (W, 0) 
-oo 

dzi / dz 2 p(b, Zl )p(b, z 2 )e izi( - k '«*- kv " ) 

(26) 



2i 



OO ./ Zl 



In the first term in Eq. f[26|) . we took into account the non-zero longitudinal momentum 
transfer, fc 7 * — k y = xbttin, see Eq. ( 12"TT) . and also neglected the t-dependence of S„ m (t) 
compated to -Fa(i). 

Using Eq. (1261) . the nuclear DVCS cross section can be expressed in the following compact 
form 

a D vcs(W,Q 2 ) = ^ J d 2 b\A{ 1 * + A-> 1 + A)\ 2 . (27) 

In order to quantify predictions of the GVMD model for nuclear DVCS cross sections, 
we introduce the ratio R cs , 

Tj _ ^DVCS^Q 2 ) , > 

where the numerator is calculated using Eq. (1271) and the complete expression for the nu- 
clear DVCS amplitude (1261) ; the denominator is calculated using only the first term (Born 
contribution) in Eq. (I26I) . 

Predictions of the GVMD model for the ratio R cs at Q 2 = 1 GeV 2 as a function of xb 
are presented in Fig. O The solid curve corresponds to 40 Ca; the dotted curve corresponds 
to 208 Pb. 

As one can see from Fig. [9], the predicted amount of nuclear shadowing at small-rcg is very 
large. Since the t-dependence of the nuclear DVCS amplitude is very steep, the dominant 
contribution to the t-integrated cross sections entering R cs comes from the t ~ t m \ n region. 
Therefore, the amount of nuclear shadowing for R cs is equal roughly twice the amount of 
nuclear shadowing for i^ pl , see Fig. [7J 

Finally, we would also like to point out that for nuclear DVCS, the ratio of the imaginary 
parts of the DVCS and DIS amplitudes, see Eq. (Tl2|) and Fig. [2j is quite similar to the free 
nucleon case. This is a mere consequence of the fact that the structure of the Q 2 -dependence 
of the ratio is essentially the same in the DVCS on the nucleon and on nuclei. 
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IV. CONCLUSIONS AND DISCUSSION 

We considered Deeply Virtual Compton Scattering (DVCS) on nucleons and nuclei in the 
framework of generalized vector meson dominance (GVMD) model. We extended the original 
GVMD model, which was applied to forward amplitudes of Deep Inelastic Scattering (DIS), 
to the non-forward t 7^ case. We introduced the W^-dependence of the DVCS amplitude 
through the Independence of the elementary vector meson-nucleon amplitude, which was 
taken to be proportional to W ' 2 at high-iy. 

We compared our predictions to the HERA data on DVCS on the proton with the fol- 
lowing results. The GVMD model describes well the dependence of the DVCS cross section 
on Q 2 , W (at Q 2 =4 GeV 2 ) and t. At Q 2 =8 GeV 2 , the ^-dependence of the cross section is 
somewhat underestimated, which can be interpreted as due to the onset of the hard regime 
beyond the soft dynamics of the GVMD model. 

We estimated the relative contribution of l/Q 2 -corrections to the DVCS amplitude and 
the DVCS cross section. We found that these corrections are large: the contribution of the 
l/Q 2 -corrections to the DVCS amplitude at t = t min is 20% at Q 2 = 2 GeV 2 , 11% at Q 2 = 4 
GeV 2 and 6% at Q 2 = 8 GeV 2 ; the contribution of the l/<5 2 -corrections to the t-integrated 
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DVCS cross section is 56% at Q 2 = 2 GeV 2 , 32% at Q 2 = 4 GeV 2 and 17% at Q 2 = 8 GeV 2 . 

We also made predictions for the DVCS amplitude and the DVCS cross section on nuclear 
targets, which are relevant for the physics program of the future Electron-Ion Collider. We 
predicted significant nuclear shadowing, which matches well predictions of the leading-twist 
nuclear shadowing in DIS on nuclei. 

Our analysis allows us to argue that the GVMD model provides a reliable parameteriza- 
tion of the DVCS amplitude and the DVCS cross section with nucleons and nuclei in a wide 
range of kinematics. At fixed values of Q 2 , which should not be too large, Q 2 < 5 GeV 2 , 
the GVMD model is applicable starting from W = 2 GeV (JLab), towards W w 80 GeV 
(HERA) and beyond (real photons at the LHC). At fixed W, the GVMD model is appli- 
cable from the photoproduction limit up to the values of Q 2 , where perturbative QCD can 
already be used, < Q 2 < 5 GeV 2 . In addition, due to the correct l/Q 2 -scaling, predictions 
of the GVMD model can be extrapolated to much higher values of Q 2 such that the range 
of applicability of the GVMD model becomes very wide, < Q 2 < 80 GeV 2 . The model 
can be applied for a wide range of t: < \t\ < 1 GeV 2 . 
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